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ABSTRACT 

We investigate the influences of the outer boundary conditions(OBCs) on the struc- 
ture of an optically thin accretion flow. We find that OBC plays an important role in 
determining the topological structure and the profiles of the density and temperature 
of the solution, therefore it should be regarded as a new parameter in the accretion disk 
CO ' model. 
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t***" ' 1. Introduction 
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, The behavior of an astrophysical accretion flow is described by a set of non-linear ordinary 

differential equations. Thus the outer boundary condition(OBC) possibly plays an important role 
in such an initial value problem. On the other hand, the complicated astrophysical environments 
make the states of the accreting gas at the outer boundary, such as its temperature and angular 

■ momentum, various. So it is necessary to investigate the role of OBC in accretion process. 

However, this problem has not received enough attention until now. In the standard thin disk 
model the differential equations are reduced to a set of linear algebraic equations that don't entail 
any boundary conditions at all. In the later works on the global solutions for slim disks the thin 
disk solution was usually adopted as the approximate OBC (Paczyhski & Bisnovatyi-Kogan 1981; 
Muchoteb & Paczyhski 1982; Matsumoto et al. 1984; Abramowicz et al. 1988; Chen & Taam 
1993). In the case of optically thin accretion, this approach does not apply any longer because 
of the importance of the energy advection throughout the disk. In this case, self-similar solution 
by Narayan & Yi(1994), or some "arbitrarily" adopted value around the self-similar solution, is 
usually chosen as the OBC(Narayan, Kato & Honma 1997; Chen, Abramowicz &: Lasota 1997; 
Nakamura et al. 1997; Manmoto, Mineshige & Kusunose 1997; Popham & Gammie 1998). They 
all concentrated on the solution under a specified boundary condition, while the role of OBC was 
not considered. Shapiro(1973) is an exception. In the context of spherically accretion, he found 
that the temperature distribution and the radiation luminosity of the disk depended sensitively on 
whether the interstellar medium was an H I or an H II region. 



In this Letter, taking the optically thin accretion as an example, we investigate the effects of 
OBC on the structure of the accreting flow. 
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2. Equations and Results 

We consider a steady axisymmetric accretion flow around a black hole of mass M{= 1OM 
throughout this Letter), employing a set of height-integrated equations widely used by the present 
workers. Paczyhski & Wiita (1980) potential is used to mimic the geometry of a Schwarzschild 
black hole. We first assume that the couple between the electrons and ions is so strong that the 
flow is a one temperature plasma. The equations read as follows: 

- A-rrrHpv = M, (1) 

v(Qr -]) = ac s r-^—, (3) 

pv dc 2 s o dp , 

c 2 s v-£= q + -q hr . (4) 



(7 — 1) dr s dr 

All the quantities have their popular meanings. The shear stress is taken to be simply proportional 

dhir 

R. 



to the pressure, i.e., shear siiess=ap d] ^f lk . We assume the flow is optically thin, so the pressure is 



dominated by gas pressure, p = jjpT. For simplicity the bremsstrahlung emission q hr is assumed 
to be the only cooling process. The viscous heating and radiation cooling rates are 

1 = ap— r— , (5) 



dlnr \ dr 

q- = 6.2 x 10 20 p 2 r 1/2 ergcm- 3 s _1 , (6) 

respectively. 

Since the accretion onto a black hole must be transonic, the solution must satisfy the sonic 
point condition at a certain sonic radius r s . Besides this, we require the flow must satisfy two outer 
boundary conditions at a certain outer boundary r out . (Note that the no torque condition required 
at the black hole horizon is automatically satisfied under our viscosity description). As M, M, a are 
given, we obtain the global solution by adjusting the eigenvalue j to ensure that the flow satisfying 
the outer boundary condition satisfies the sonic point condition as well. 

At r ou t, the two outer boundary conditions we imposed are T — T ou t and A(= -7- = 



'p/p 

A out . This set of boundary conditions are equivalent to (T out , ri out ) according to eq. (3) because j 
is not a free parameter. We assign (T out , A ou t) to different sets of values and investigate their effects 
on the global solutions. The results are as follows. 

As M, a, and r out are given, we find that only when T out and A out are within a certain range 
do the global solutions exist. For example, for the global solutions with M = 10~ 3 Me ct = 10~ 2 
and r out = 10 3 r g (Here the Eddington accretion rate reads Me = 4 ^f M and r g is the Schwarzschild 
radius), the values of T out are close to the virial temperature T v i r i a i(= (7 — l)GMm/Kr out ), T out ~ 
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0.1 — lT v i r i a i. This range changes slightly when we modify the values of M, a and r out , and when 
we consider different radiation mechanisms. 

The solutions are remarkably different under different OBCs. In terms of the values of T out , 
the whole parameter space spanned by T ou t and A ou t can be divided into two regions. The solutions 
located in the low-temperature region(referred to as type I) have the following features, as the solid 
line in Fig. 1 shows: 1) the maximum feasible value of A ou t is close to 1; 2) a maximum exists in 
the solution's angular momentum (/ = ilr 2 ) profile; 3) the sonic radius r s is always small; and 4) 
the Bernoulli parameter Be, which is equivalent to \v 2 + ^zjc^ + — ^ s negative at large 
radii. The solutions located in the high-temperature region can be further divided into two types 
in terms of the values of A ou t- When A ou t is smaller than a critical value A ou t,crit (referred to as 
type II), the sonic radii of the corresponding solutions are still small, but the values of Be become 
positive. Besides these, the angular momentum and the surface density(S = 2pH) profiles and 
the topological structures (Mach number profile) differ greatly from those of type I, as shown by 
the dot-dashed line in Fig. 1. When A is greater than A cr i t (referred to as type III), r s suddenly 
becomes very large, as Table 1 and the dashed line in Fig. 1 show. This is a new kind of solution, 
since all previous works on ADAFs global solutions only obtained solutions with small values of r s . 
This kind of solution is difficult to be found because the parameter space where it is located is too 
narrow. 

It is interesting to note that such transition of the values of r s is similar to that in the adiabatic 
case(Abramowicz & Zurek 1981). In that case, Be and I are two constants of motion. When Be 
lies in a certain range, r s will jump from a small value to a large one when I decreases across a 
critical value l c , and the corresponding accretion patterns are called disk-like(small r s , large I) and 
Bondi-like (large r s , small I) ones respectively. Our result indicates that such transition still exist 
although the flow becomes viscous (note that the increase of A corresponds to the decrease of the 
angular momentum). 

In general, modifying A ou t(i.e., modifying O, ont ) has much smaller effects on the solution than 
modifying T out does. The latter has a great effect on the topological structure and the profiles 
of the angular momentum and surface density of the solutions. But its effect on the temperature 
profile rapidly weakens with the decreasing radii, as the temperature profile plot in Fig. 1 shows. 
We'll discuss this problem later. 

The above calculation indicates that, besides M and a, OBC is also a crucial factor to deter- 
mine the structure of the accretion flow. 

Present theoretical and observational works seem to support the two temperature plasma 
assumption with the ions being much hotter than the electrons. In the following we consider 
whether the effects of OBC are still important in this case. Comptonization of the bremsstrahlung 
photon is also included in the cooling processes. The pressure and the energy equations should be 
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replaced by, 

pkTi pkT e 

P = Pi+Pc = 1 , (7) 

pv (^p- + pid (-\ \ =q+ - q ie , (8) 



3 m i 1 )' +(0e + 0i) 1/2 
Qie = -—n eni o T c\xik (m - kT e ) — , (11) 



V dr VP 

PU + Ped (J) ) = qie ~ %c ' (9) 
where \x x and /U e denote the ion and electron molecular weights respectively, Ei,e e are the internal 
energies of the ion and electron per unit mass of the gas. Since for the temperature of interest, the 
ion is always unrelativistic while the electron is transrelativistic, following Shapiro(1973), we adopt 
the following: 

3 kT{ kT e 

Si = , £ e = 2 . (10) 

2 m m ti Vemn 

qi e in eq. (9) is the energy-transfer rate from the ions to electrons per unit volume duo to Coulomb 
collisions(Dermer, Liang & Canfield 1991): 

where InA = 20 is the Coulomb logarithm, and 6 = kT/mc 2 is the dimensionless temperature. 
The thermal bremsstrahlung amplified by Comptonization q^ c reads( Rybicki & Lightman 1979) as 
follows: 

%c = CrA (12) 

A = m a x{^\n 2 [x ff ^ln(xj} )],1}, (13) 
x f f = 2.8 x W 12 T~ 7 / 4 pH 1 / 2 . (14) 

Eqs. (7)-(14) together with eqs. (l)-(3) constitute the set of equations describing a two 
temperature accretion flow. In this case we must supply three physical quantities at r out , e.g. 
T out ,i,T outtC and A out (= vj^fpfp). Obviously, for a given T out)i , the upper limit of T outiC is T out) i 
while its lower limit is determined by q- lc = q hc , i.e. all the energy transferred to the electrons is 
radiated away. 

We solve the above set of equations under different OBCs and find that the results are quali- 
tatively the same, as shown by Fig. 2. In this case we still find three types of solutions, the effect 
of A ou t is still unimportant for an individual type of solution, and the surface density profiles still 
differ greatly for the solutions with different OBCs. Compared with the corresponding solution in 
Fig. 1, the value of Be of the type I solution (denoted by the solid line in Fig. 2) decreases greatly. 
This is due to the inclusion of Comptonization together with the high surface density of the flow. 



A remarkable new feature is that the electron temperature profiles of the solutions with different 
OBCs(primarily T outje ) strongly differ from each other throughout the disk, while for the ion the 
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discrepancies rapidly lessen with the decreasing radii from r out like a one temperature plasma. 
This is because the electrons are essentially adiabatic in the present low M case, i.e., both g; e and 
qT are very small compared with the two terms on the left hand side of eq. (9), so the electron 
temperature is globally determined therefore the effect of OBC on T e persists as the radius decreases. 
But for the ion, as well as the one temperature plasma, the local viscous dissipation in the energy 
equation plays an important role thus the energy equation for ions is much more "local" than that 
for electrons, so the effect of OBC weakens rapidly with the decreasing radii. Due to the same 
reason, we predict that OBC will be less important in determining the temperature profile of an 
optically thick accretion flow, since both the viscous dissipation and the radiation losses terms in 
the energy equation are important thus the temperature will be locally determined. But the flow 
should still present OBC-dependent dynamical behavior in , for example, the surface density and 
angular momenta profiles and the topological structure. 

3. Conclusions 

In this Letter, taking the optically thin one temperature and two temperature plasma as 
examples, we investigate the role of the outer boundary condition in the accretion process and find 
that it plays an important role. Our results indicate that in either case the whole solutions can 
be divided into three types in terms of their OBCs. The topological structures and the profiles 
of the angular momenta and the density are strongly different among the three types. As to 
the temperature we find that for a two temperature plasma the discrepancies of T outiC among the 
solutions with different OBCs will persist throughout the disk, while for ions and a one temperature 
plasma the discrepancies of the temperature rapidly lessen away from the outer boundary. Such 
discrepancies among the dynamic quantities will further induce the discrepancies of the radiation 
luminosity and the emission spectra of the accretion disc. Therefore we argue that OBC should be 
regarded as a new parameter in the accretion model. 

F.Y. thanks Dr. Kenji E. Nakamura for helpful discussions and the anonymous referee for 
useful suggestions to improve the presentation of this paper. This work is supported in part by the 
National Natural Science Foundation of China under grants 19873005 and 19873007. 
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Table 1. 

Parameters, OBCs and the sonic radii of some solutions 
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Fig. 1. — Solutions for one temperature global solutions with different OBCs for M = 10 3 Me 
and a = 10 -2 . The solid, dot-dashed and dashed lines represent (T out , A out )=(2 x 10 8 K, 0.4), (3.6 x 
10 9 K,0.08) and (3.6 x 10 9 K, 0.107) respectively. The units of £,T are g cm' 2 and K, r, Be and I 
are inc = G = M = l units. Mach number is simply defined as v/c s . The upper left-hand plot 
represents the ratio of the advected energy to the viscous dissipated energy. 



-9- 



CO 




1 2 3 1 2 

log(r) log(r) 



Fig. 2. — Solutions for two temperature accretion flows with different OBCs for M = 10 3 Me 
and a = 10~ 2 . The solid(type I), dotted(type II), dashed(type III) and dot-dashed (type II) lines 
represent (T outji , T outjC , A)=(6 x 10 8 K, 5 x 10 8 K, 0.8), (3.2 x 10 9 K, 2.62 x 10 8 K , 0.1), (3.2 x 10 9 K, 3 x 
10 9 K,0.1) and (3.2 x 10 9 K,3 x 10 9 K, 0.088), respectively. The units are same as those in Fig. 1. 
In the dotted line T outjC is determined by setting q m = q^ c at r out . 



